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Abstract 

Stock price change in financial market occurs through transactions 
in analogy with diffusion in stochastic physical systems. The analysis 
of price changes in real markets shows that long-range correlations of 
price fluctuations largely depend on the number of transactions. We 
introduce the multiplicative stochastic model of time interval between 
trades and analyze spectral density and correlations of the number 
of transactions. The model reproduces spectral properties of the 
real markets and explains the mechanism of power law distribution 
of trading activity. Our study provides an evidence that statistical 
properties of financial markets are enclosed in the statistics of the 
time interval between trades. Multiplicative stochastic diffusion may 
serve as a consistent model for this statistics. 

Keywords: statistical analysis, financial markets, stochastic mod- 
elling. 

1 Introduction 

Complex collective phenomena usually are responsible for power-laws which 
are universal and independent of the microscopic details of the phenomenon. 
Examples in physics are numerous. Power-laws are intrinsic features of the 
economic and financial data as well. The aim of this contribution is to 
analyze a relation between the origin of the power law distribution and the 
power-law correlations in financial time series. There are numerous studies 
of power-law probability distributions in various economic systems fil ||, 
[|, ||]. The time-correlations in financial time series are studied extensively 
as well [H, ||, The random multiplicative process built into the model 
of wealth distribution yields Pareto style power law Q. The generalized 
Lotka Volterra dynamics developed by S.Solomon and P.Richmond is in 
the use for various systems including financial markets |l] . However, 
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these models generically lead to the non-universal exponents and do not 
explain the power-law correlations in financial time series ||. Recently we 
adopted the model of 1/f noise based on the Brownian motion of time 
interval between subsequent pulses, proposed by Kaulakys and Meskauskas 



|j, 1C, 11, [T^] to model share volume traded in financial markets ]L_: 
. The idea to transfer long time correlations into stochastic process of 
the time interval between trades or time series of trading activity is in 
consistence with the detailed studies of the empirical financial data 
and fruitfuly reproduces spectral properties of financial time series 
However, the investigation of the model revealed, that the simple additive 
Brownian model of time interval between trades failed to reproduce power- 
law probability density distribution (pdf) of trading activity. On the other 
hand, several authors showed empirically that the fluctuations of various 
financial time series possess multifractal statistics [14, 15, 1€] . Therefore 
we introduce the stochastic multiplicative model of time interval between 
trades and analyze the statistical properties of the model trading activity 
numerically and analytically. We show the consistence of the approach with 
the results of statistical analyzes of the empirical financial time series. This 
sheds light on the relation between the power-law probability distribution 
and the power-law correlations in financial time series. 



2 Multiplicative time interval model 

Following our previous model |[3| we consider a signal I(t) as a sequence 
of the random correlated pulses 

k 

where is a contribution to the signal of one pulse at the time moment 
tk, for example, a contribution of one transaction to financial data. Signal 
(1) was introduced by Kaulakys and Meskauskas for the modelling of 1// 
noise and could be used in a large variety of systems with the flow of point 
objects or subsequent actions. In the statistical data analysis we usually 
deal with the discrete time series in equal time intervals Integrating 
the signal I(t) in the subsequent time intervals of length we get discrete 
time series and, in analogy with financial time series, we call it the volume 
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tj+T d 



Vj= / I(t)dt= Yl 9k 



tj = JT d . 



(2) 



tj <t k <tj+T d 



We can define the number of trades Nj in the time interval tj -r- (tj + Td) by 
the same formula (2), when = 1 , Nj = Vj . The power spectral density 
S(f s ) of the signal Vj is defined as 



where f s = (s — 1)/T , T = r^n. In this paper we investigate statistical 
properties of the time series Nj , when the sequence of event time if. is gen- 
erated by the multiplicative stochastic process. There are some arguments 
for this choice of interevent time model. First of all the multiplicativity is 
an essential feature of the processes in economics [|, |, H H . Multiplicative 
stochastic processes yield multifractal intermittency and are able to produce 
power-law pdf. Pure multiplicative processes are not stationary and artifi- 
cial diffusion restriction mechanisms have to be introduced. Let us start our 
analysis introducing the multiplicative process with the logaritmic diffusion 
restriction for the large deviations of from the mean value = 1 

tk+l = tk + Tk+l 



Here the interevent time fluctuates due to the external random perturba- 
tion by a sequence of uncorrelated normally distributed random variables 
{£k} with zero expectation and unit variance, where a denotes the standard 
deviation of the white noise and 7 << 1 is a damping constant. The choice 
to introduce multiplicative logaritmic diffusion damping 7^. In 7^ retains 
symmetry in logaritmic scale and multiplicative nature of the process for 
small values of 7^7^ . Long time stationary pdf of r& is a Lognormal 
distribution 




(3) 




(4) 




(5) 
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We derive the distribution (5) from the solution of the generalized Langevine 



equation 17] 



x = F(x) + G(x)r](t), 
(V(x) 



p(x) = C exp 



\° 2 



dV _ F 
~dx~ ~ CP' 



+ lnG(x) 



(6) 



where {rj(t)} = and (r/(t)r](t')} = a 2 5(t - t'). Note that Lo gnormal 
distribution is invariant to the change of variables: n = 1/r and in the wide 
range of argument values |lnr| << cr/7 1 / 2 is equivalent to the power-law 
~ 1/r. Nevertheless, we have to point out that lognormal distribution of r k 
is not appropriate to reproduce the power-law distributions of variables in 
financial markets. We expect that iterative stochastic model (4) producing 
universal lognormal distribution will be useful for other applications. We 
base our model on the generic multiplicative process 

r k+1 = r k + jt^' 1 + r^ae k , (7) 

as the most appropriate seeking to reproduce power-law probability density- 
distributions of variables in financial markets. In this paper we will use a 
very simple r k diffusion damping model assuming that < 77 < 1 and < 
7 << a. More sophisticated damping models can be introduced in further 
developments of this approach. The most natural choice is to assume fi = l, 
giving pure multiplicativity of 77. We decided to keep this parameter in 
our notation seeking to sustain generic definition. In the k scale Eq. (7) 
defines generalized Langevine equation (6) with F(t) = 77"^ 1 and G(r) = 
Tu ■ Note that the solution (6) of the generalized Langevine equation is 
based on the Stratanovich convention. Our numerical analysis confirms the 
preference of the Stratanovich convention in comparison with Ito convention 
regarding discrete iterative equation (7). The stationary pdf of r k defined by 
iterative relation (7) can be derived from the solution (6) of the generalized 
Langevine equation, 

W = (8) 

where the normalization constant Crr) can be defined from the integral 
1 

J Pk(r)dT = 1. The stationary distribution function for the number of 


trades N k in the time interval can be derived from (8). Simple change 
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of variables N = t^/t in (8) having in mind that Pk{r)dT = Pk(N)dN , 
gives probability distribution function of N 



where the time interval is included in normalization constant Cn- Note 
that variance (N?) — (N^) 2 diverges when 2j/a 2 < /j, + 1. We will see later 
that the transition from the stationary regime to the non-stationary one is 
essential for the appearance of 1/f power spectral density. Distribution 
function (9) can be transformed to the time scale taking into account that 
every step in the iteration equation (7) is equal to the time interval r^. In 
the real time scale the pdf of N takes the form 

Here CV t is a new normalization constant. 



3 Numerical results and discussion 

We have already introduced the multiplicative stochastic process as a model 
for the time interval between trades in financial markets. This yields the 
power-law distribution for the number of trades per constant time interval 
Td, equation (10). The exponent of cumulative power-law distribution a 
for pure multiplicative model, \i = 1, depends only on the parameter v = 
2 7 /a 2 , i.e., 

a = l + 2-f/a 2 . (11) 

Parameter v defines the ratio of damping constant 7 to the coefficient of 
stochastic diffusion \o 2 . In this chapter we will analyze spectral properties 
of the model and related autocorrelation of trading activity N. The model 
presented can be easily calculated numerically using equations (1), (2), (3) 
and (7). We smooth the power spectral density with standard moving 
average procedure. We present the numerically calculated power spectral 
density S(f) of N for v = 2.4; 2.0 and 1.5 on the Fig. 1. (a), (b) and (c) 
respectively with a = 0.1 and = 10. For the values of parameter v > 2, 
we observe the power spectral density S(f) ~ 1/f 1 ' 5 . For v ~ 2 we obtain 
S(f) ~ 1// and we get the slope of S(f) ~ 1/f - 5 when v < 2. 
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Figure 1: Power spectral density of trading activity N (N PSD) versus 
frequency / calculated from the model described by Eqs. (1), (2), (3) and 
(7) with fj, = 1. The sinuous curves represent results of the numerical 
simulations smoothed with standard moving average procedure over 43 
points with a = 0.1 and Td = 10, while the straight lines represent power- 
law power spectral density S(f) ~ 1//^ with the various slopes [3. For 
7 = 0.012, = 1.5 (a); for 7 = 0.01, = 1 (b) and for 7 = 0.0075, (3 = 0.5 
(c). 
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Figure 2: The same as in fig. 1 but with parameters a = 0.032 and Td = 
100 or 200. The straight lines represent power-law power spectral density 
S(f) ~ with the various slopes (3. For 7 = 0.0012, (3 = 1.5 (a); for 

7 = 0.001, (3 = 1 (b); for 7 = 0.0008, (3 = 0.5 (c). 

Numerical results confirm that 1/f noise occurs when stochastic process 
defining the trading activity experience the transfer from the stationary 
to the non-stationary regime. This general rule is probably applicable to 
the various models of the time intervals between pulses and is in agreement 
with the theory of self-affine fractals, see |l8| . It should be noted that 
transition from slope (3 = 1.5 to (3 = 0.5 of the power spectral density is 
very sharp with the change of parameters v or 7 less then twice. Calcula- 
tions with different values of parameters confirm a universal nature of the 
relation between the slope of the power spectral density (3 and parameter of 
multiplicative diffusion v . This relation is independent of in wide range 
of values. We present our numerical results of power spectral density for 
different sets of parameters in Fig. 2. The slopes (3 = 1.5; 1.0an<i0.5, 
of power spectral density are related with the corresponding values of 
parameter v = 2.4; 2.0andl.6, as in Fig. 1. Numerous calculations confirm 
that the slope of power spectral density f3 is defined by z^, exhibiting sharp 
transition from (3 = 1.5 to (3 = 0.5, when v ~ 2.0. This transition is related 
with divergence of iV/% variance (N?) — (N^) 2 defined from (9). 



In the previous section we derived the power-law probability distribution 
function of N, Eq. (10). Let us compare the exponent of cumulative 
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Figure 3: The exponent a of power-law distribution function for N versus 
multiplicative diffusion parameter v calculated numerically from the model 
described by Eqs. (1), (2), (3) and (7). Points represent various realizations 
with the parameters a = 0.05 4 0.1, 7 = 0.001 4 0.015, r d = 10 4- 100. 
Straight line represents the least square approximation to the numerical 
data a = 0.96 + 0.97u. 

distribution from (10), a = 1 + v , with numerical results in the interval of 
v values 1.5 4- 3.5 related with the slope (3 change. In Fig. 3. we present 
numerically calculated a as a function of v for various values of a, 7, and r^. 
The least square linear approximation of numerical results, straight line, 
yields that a = 0.96 + 0.97u. It shows that the probability distribution 
function (10) derived from the generalized Langevine equation (6) describes 
the power-law distribution of N very well and serves as an evidence that 
the Stratanovich convention is appropriate for the multiplicative stochastic 
model (7). For v ~ 2, corresponding to (3 = 1., we get the value of a ~ 3. 
This is in good agreement with empirical data of power-law distribution of 
trading activity N in the financial markets a ~ 3.4 ||. 
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Numerical results confirm that multiplicative stochastic model of the time 
intervals between trades in the financial markets is able to reproduce the 
main statistical properties of trading activity N. First of all, even in a very 
simple multiplicative diffusion restriction < Tk < 1, the model reproduces 
the power-law distribution function with a ~ 3. The power-law exponent 
is related with the parameter v = which defines the distribution 

functions of time interval between trades r, and of trading activity N as 
well as the slope (3 of its' power spectral density S(f) ~ 1// /3 . This sheds 
light on the origin of universal nature of the power-law distributions of such 
variables as trading activity in financial markets N and its' power spectral 
density S(f). 

4 Conclusions 

We have introduced multiplicative stochastic model of the time interval r 
between trades in the financial markets Eq. (7) with the diffusion restriction 
in the interval < Tfc < 1. This defines stochastic fluctuations of trading 
activity N, Eq. (2), and reproduces its' power-law distribution of probabil- 
ity and the slope of power spectral density S(f) ~ I//' 3 - The comparison 
of the empirical data from the financial market analysis with the model 
proposed implies that markets function in the area of parameters, where 
convergence of variance disappears and 1/f noise of N occurs. This 
rather a sharp transition to the non-stationary regime of fluctuations 
probably is the indispensable feature of the financial markets. This relates 
the exponent of the power-law distribution to the slope (3 of power spectral 
density into universal interdependence. It implies that the main statistical 
properties of the financial markets have to be defined by the fluctuations of 
time interval between the trades accumulating the power-laws and the long 
time correlations. Further empirical analysis of r statistics and the adjusted 
specification of the model is desirable. We do expect that multiplicative 
model of time interval between the trades with more specific diffusion 
restriction conditions and more precisely adjusted fx lies in the background 
of the financial market statistics and can be very useful in financial time 
series analysis. 
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